Thermal fluctuations of bouncing cosmology revisited 
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In this note, we revisit the thermal fluctuations generated during bouncing cosmology, taking 
Unruh effect into account. We find that due to the additional effect on temperature, the dependence 
of power spectrum on k will get corrected with an indication of blue tilt at large k region, which is 
in consistent with the case of vacuum initial conditions. 
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INTRODUCTION 

The well-known observational data of Cosmic Mi- 
crowave Background anisotropy as well as Large Scale 
Structure [l| shows that the structure of our universe, 
such as galaxies, temperature fluctuations and so on, 
comes from the primordial perturbation at early times 
with a nearly scale invariant power spectrum, as success- 
fully predicted by Standard Big Bang Theory I n 
common sense, this perturbation was initiated from vac- 
uum quantum fluctuations, which was then stretched out 
of the horizon. It was noticed, however, that the vacuum 
initial condition is by no means the only way to gener- 
ate cosmic perturbations and form structure. Thermal 
fluctuation, of which fluctuations comes from a thermal 
ensemble, has been argued to be the alternative origin to 
source the cosmic structure 0. 

Unfortunately, as was concluded in [7] and summa- 
rized in the following paper Q, that the scale invariant 
power spectrum can not be generated in normal infla- 
tionary stage. Due to this reason, people begin to find 
new physics in order to save the scale invariance, Among 
which are string gas cosmology [9, LL0( , non-commutative 
inflation scenario [11 1, h olographic cosmolog y |l 2|. loop 
quantum cosmology J A and so on. See also [14H171 ]. 

Recently, people paid more and more attention to the 
scenario of non-singular bouncing cosmology, which has 
the advantage of getting rid of the initial singularity 
which plagued the Standard Big B ang or inflation the- 
ory of cosmology ([HHHj], also see |32| for a review and 
the references therein). Starting from a large and cold 
initial state, a bounce scenario can also lead to scale in- 
variant power spectrum of perturbation, It was then ar- 
gued in [33| that thermal fluctuations is not only possible 
to occur, but also able to obtain scale invariant power 
spectrum as the case of vacuum initial conditions. The 
fluctuations can either be ensemble of particle gas, or 
holographic radiation, or even string gas. Imprints of 
non-Gaussianities could also been obtained due to ther- 
mal fluctuations. 

In this paper, we revisit the thermal fluctuations gen- 
erated during bouncing cosmologies, taking for exam- 



ple that the fluctuations as normal particle radiation of 
which the equation of state w r = 4. When the fluctua- 
tion moves with varying velocity, i.e. with an accelera- 
tion (which is a very common case) , then an Unruh effect 
should be taken into account. It will cause an additional 
"Unruh" temperature which may give a correction to the 
original temperature T r . We have calculated that it will 
scale as T r times some power law of k. For the origi- 
nal case which can get scale invariant power spectrum, 
the power law index is positive and will give a enhance- 
ment at large k region (small scales). It will cause an 
indication of blue tilt in the power spectrum, which is 
in consistent with what happens in the case of vacuum 
initial conditions 31 1. 

The paper will be organized as follows: we first remind 
the reader the main steps in calculating thermal fluctua- 
tions in bounce cosmology in the next section, and then 
give our results on correction to the spectrum index by 
taking the Unruh effect into account. The last paragraph 
is the summary. We leave the basic conjecture for back- 
ground evolution and the formulae for perturbation in 
the Appendix A. 



THERMAL FLUCTUATIONS IN BOUNCING 
COSMOLOGY 

In this section, we review the thermal fluctuations in 
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bouncing cosmology, as is already mentioned in [31 
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In order to obtain the spectrum at Hubble crossing at the 
expanding phase, we need to calculate it first in the con- 
tracting phase, then transform it into expanding phase 
using matching conditions at the bouncing point. We put 
the general argument of the metric perturbation evolu- 
tion and matching conditions to Appendix A. First of all, 
from a perturbed Friedmann-Robertson- Walker (FRW) 
metric, 

ds 2 = a(ry) 2 [(l + 2$)d7y 2 - (1 - 2$)d^ 2 ] , (1) 

we write down the first perturbed Einstein Equation as: 

-3-H($' + -H$) + V 2 $ = 47rGa 2 (5p. (2) 

In the above two formulae, a(rj) is the scale factor in 
terms of the conformal time rj = J ^4|j, $ is the Newto- 
nian potential, H = is the Hubble parameter in the 
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conformal time, 8p is the perturbation of energy density, 
and prime denotes for derivative with respect to r\. As- 
sumption is taken that we are in conformal Newtonian 
gauge and there is no anisotropic stress 34(. At Hubble 
radius crossing all the three terms on the left-hand side of 
the above equation are of the same order of magnitude, 
so modulo a constant of the order one, we can write down 
the Fourier space correlation function of $ as: 



|$(fc)| 2 = 167r 2 G 2 fc- 4 a 4 < Sp(k) 2 > 



(3) 



In the system of thermal equilibrium, the energy 
density fluctuations are determined by thermodynam- 
ics. The correlation function of Sp in momentum space 
are correlated to that in position space by the relation 
< Sp(k) 2 >~ fc~ 3 < 5p(x) 2 >, where the latter can be 
represented as: 



< 5p{xf >= C V {R) 



(4) 



where Cy (R) = d< ^p > is the heat capacity in a sphere of 
radius R. 

We assume that the fluctuations consist of an ensem- 
ble of hot particle gas with constant equation of state 
w r . Then from the continuation function and the second 
thermodynamics law, one can derive the energy density 
and temperature which scale as: 



Pr ~ 

respectively. This gives 



-3(l+«v) 



-3w r 



Pr 



T, 



1+^7- 



(5) 
(6) 



(7) 



and thus 



, r l \2 v. -^r dp r 

<Sp{x) > = mW r 

q 2-1 — *- 

~ R T r Wr 

^ £—3(l+p) — 6pw r 



(8) 



where in the last formulation we have used R = -tj . Com- 
bining Eqs. © and (JSJ, one can obtain the dependence 
on t and fc of |$| 2 that is: 



1$ 



|2 ^ fc 3 < 5p{x) 2 > ^ -3.1-3p-6pu 



and 



H 



$ ~ k 5 f 2 



(9) 



(10) 



At the Hubble crossing, where fc* = , we have 

l-3p(l + 2 mr ) 

^ , (11) 



thus the power spectrum at Hubble crossing time will 
be (See Appendix A for more details): 

L.3 l-3p(l + 2tu r ) 

= £-l$l 2 ~ &; ^ . (12) 

Z7T 

From above we can see that the condition for the scale 
invariance of the spectrum is 



1 — 'ip — 6pw r = 



(13) 



and for the case that the fluctuation behaves like normal 
radiation, w r = \, it is straightforwardly obtained that 

1 • 7 

p = g and the equation of state of the background w = ^ . 



CORRECTION FROM UNRUH EFFECT 

The Unruh effect, which was first described by S. 
Fulling in 1973, P. Davies in 1975 and B. Unruh in 1976 
(35^ predicts that an accelerating observer will observe 
non-zero temperature where an initial observer would sec 
none. In other words, from the viewpoint of the observer, 
particles which is accelerating will look like a warm gas 
which is in thermal equilibrium with non-zero tempera- 
ture [3^ . As a consequence, it will cause the black body 
radiation as Unruh radiation, which can be viewed as the 
near-horizon form of Hawking radiation (see e. g. (37j 
for a comprehensive review and references therein). 

We consider the Unruh effect by taking the velocity 
part of the fluctuation into account. The — i component 
of the perturbed Einstein Equations is: 



+ = A-KGa{p + P)5ui , 



(14) 



where ", i" denotes the covariant derivative with respect 
to the 3-dimensional metric 7^ , and Sui is the velocity 
of perturbation in i direction. Ignoring the anisotropy of 
the perturbation and assuming Sui — 8u2 — Su-j = 8u, 
and presenting the perturbation in its Fourier form, we 
have: 



|*u(fc)| 



I 4nGa(p + P) 



(15) 



It is useful to define the comoving curvature perturba- 
tion: 



w-n 



7 ($' + h$) 



which satisfies the equation 



C" + 2-c + c 2 fc 2 c = o, 



(16) 



(17) 



where z = —j=p==. Furthermore, by assuming that in 
the whole process the non-adiabatic perturbation can be 
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ignored 1 , we have: 

C' = 



fc 2 ($' + H$) 
l2irGa 2 {p + P) ' 



(18) 



which shows that ( is conserved at large scales. 

From (TT5|) and (|18[) . and considering the velocity per- 
turbation in position space, Su(x) ~ ki5u(k), we have 



|5u(x)| = 3afc5|C|' 



(19) 



which relates the dynamics of the fluctuations to the met- 
ric perturbation. In general, the fluctuation is not mov- 
ing uniformly and with an acceleration a u . Thus from 
above, we can define the Unruh temperature of the per- 
turbation to be: 



2tt 



2ir 



1 

2^ 



(3a 2 ^|C|)\ (20) 



where a dot means derivative with respect to cosmic time 
t and we used 4- = -4- . 

at aarj 

Since all the terms in (1161) are of the same order of 



magnitude, then up to a coefficient of order one, we have: 
T u = i(3a 2 fci|6|)-, 

Z7T 

~ ^kHk-h 1 -^^)) , (21) 



where Ea. (ITU|) is taken into account. Then we have: 

T u ~ k-H*Tr- apu ' r ~ k-H*TrT r . (22) 



At the Hubble crossing, we replace t by kp- 1 to get: 



(23) 



This is the relation between the temperature of the ther- 
mal gas fluctuation and its Unruh correction, and the re- 
lationship depends on the equation of state of the back- 
ground 2 . For example, for the w = | case of which a 
scale invariant power spectrum can be obtained, we have 
2 - i and thus T u ~ k*T r . 



P 



3(1+™) 



This is our main result in this paper. We conclude that 
the Unruh Temperature caused by thermal fluctuations 
scales as the original temperature times a term that has 
positive power law of k. This means that in the region 
of large k, the Unruh effect will get larger and have a 
magnificant backreaction on the fluctuation. 

Next We will roughly estimate how the backreaction 
would be. From the above relation (l23l) between T r and 



T u , we can similarly parameterize T u 



where one 



can calculate easily that at the Hubble crossing w' r — 
w r — At the large k region where T u dominates, one 
can get the dependence of the power spectrum on k as 

1 — 3p — Gpw r 

what has be done at the above, V${k) ~ k ~ . For 
our case p = | and w r = 5, we can get w' r = |, and 

thus V^ik) ~ hi. From this we can see, at the large k 
region, the power spectrum will get an indication of blue 
tilt, which is in consistent with the case of vacuum initial 
conditions. 



SUMMARY 

In this short note, we reinvestigated the thermal fluctu- 
ations generated during bouncing cosmology, taking into 
account the corrections from Unruh effect. Our calcu- 
lation shows that for the original case of scale invariant 
power spectrum, the additional correction leads to a blue 
tilt at large k region, which is in consistent with the case 
of vacuum initial condition as previously discussed. 

Unruh effect is of great importance in its own right 
and has been applied as a tool to investigate other phe- 
nomena such as the thermal emission of particles from 
black holes 38 , 3(| and cosmological horizons ■ How- 
ever, its applications on thermodynamics of cosmological 
perturbations is not clear yet. Whether it will give rise 
to more corrections or modifications to the early uni- 
verse is interesting and worth studying. It can also been 
applied to the Holographic Cosmology, cf. the recently 
released paper by E. Verlinde [H| who argues that the 
Einstein Gravity is emergent from thermodynamics and 
holographic principle. These extensions of this project 
are expected to be done in the future work. 



It is a reasonable assumption and can be realized easily in bounce 
cosmology if we require the bounce occurs fast enough. For ex- 
ample, in two fields bouncing model, the short duration of the 
bounce relates to small ratio between masses of light and heavy 
fields, which can then suppress the isocurvature perturbation. 
This can be seen in Ref. |3ll . 

One may be worry about whether the amplitude of Unruh tem- 
perature can be large enough to really affect the original tem- 
perature T r . We point out that the analysis in this paper is 
quite qualitative, and the specific amplitude is determined by 
the initial condition of background evolution of the model. Given 
proper initial conditions, it is not difficult to let the corrections 
comparable to T r . We thank Yi Wang for reminding us this 
point. 
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Appendix A 

The appendix aims at clarifying the basic founda- 
tions of background and perturbation evolution for our 
analysis. In the framework of Friedmann-Robertson- 
Walker(FRW) geometry and taking the assumption of 
some constant equation of state w for the background, 
the scale factor of the universe scales as: 



where 



P 



3(1 + w) 



and the Hubble parameter becomes: 
H 



V 
t 



(24) 



(25) 



(26) 



Thus from the conditions for the comoving wave mode k 
at Hubble crossing, we have 



K = a(U)H(U) ~ t*~ 



(27) 



Then we briefly review the perturbation evolution 
through the bounce. From the perturbed Einstein equa- 
tion, we can get the equation for $ as: 



2H 3 - mw 



w-w 



U" x . UU" - 2H' 2 x 

$ H — — TT^$ 



H 2 -H' 



= . 

(28) 

The equation has a general solution of $(&, rj) which can 
be parameterized as: 



$(/c,?7) ± = D(k)± + S(k)±fj 



-2;/ 



(29) 



where v = 2 (i+3w) an< ^ re f ers to modes in the expand- 
ing (contracting) phase. Applying the Hwang- Vishniac 
[42( | (Dueruelle-Mukhanov [43() conditions which match 
modes before and after the bounce, we can get the dom- 
inant mode for the expanding phase turns out to be 44- 



D+ = 0(1)D. 



0(l)(f) 2 S_ 



(30) 



The spectrum in the expanding phase can only be given 
by the £)_ mode [33j |. which indicates that 



(31) 
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